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, Five einonwtors of th<> linen thoniMiviM-ous Hew ovor tho lionios-eneous l»;n-knn>uinl derived for the 

quasi-plane geometry of the flow. The corresponding projectors are calculated and applied to get the 
nonlinear evolution equations for the interacting vortical and acoustic modes. Equation on streaming cased 
by arbitrary acoustic wave is specified. The correspondence to the known results on streaming cased by 
■ quasi-periodic source is traced. The radiation acoustic force is calculated for the mono-polar source. 

1 INTRODUCTION 

Traditionally, streaming is referred as a localized mean flow caused by an acoustic source. This flow is created 
by the Reynolds forces, non-zero time-averaged values of the quadratic acoustic terms arising in the equation 
fH , of conservation of momentum. There exist extensive reviews on this subject [1-3]. 

Steady acoustic streaming is described by the mean Eulerian velocity ( v ) , or the mean mass flow velocity U , 

where U is defined by U — ^ - — V + ( Pa ^ a ) Here, V is a slow component of velocity, namely the Eulerian 
velocity of the streaming, index a labels acoustic wave, po is undisturbed density. Square brackets designate 
averaging over time whose interval is much shorter than the transient time of streaming and much longer than 
a period of a sound wave. It is pointed, that the difference between U and V is roughly proportional to the 
acoustic intensity vector and is small compared with any of mean velocities [4] . From the continuity equation 



it follows V [/ = if we assume ( §£ > = 0. 



The recent discussion of acoustic streaming [5,6] suggest that there is important unresolved issue concerning 
acoustic streaming, namely the effect of fluid incompressibility. The inconsistency of this point is that fluid is 
supposed to be incompressible while the acoustic wave casing streaming may propagate only over compressible 
medium. The next, the conservation of energy usually is not considered. It is reasonable since excluding of 
one variable (density) reduces an initial system. It has been proved that the effect of the heat conductivity 
could not be discarded in a study of temperature variation associated with the streaming [6] . Actually, this 
approximation is well understood in a typical liquid like water but should be revised for other liquids [7] . That 
is a reason to start from the full system of conservation laws including the energy balance equation and the 
continuity equation. 

The procedure of the temporal averaging needs a continuous acoustic wave of many periods as a source, 
thus the well-known results on acoustic streaming apply only for this type of acoustic source though quite 
realistic. There are clear inconsistency concerning the very procedure of the temporal averaging. As it has been 
mentioned before, an equality Vf7 = is correct only if ^§f^> = . Actually, the overall field includes not 

only acoustic and vortical modes, but also of the entropy mode. The entropy mode caused by acoustic wave is 
known as a slow decrease of density and isobaric increase of temperature of the background of acoustic wave 
propagation. Recently, the importance of change in density of surrounding (not only increase in temperature) 
was underlined in the review [8]. If nonlinear interactions are accounted in a proper way, for density of the 
entropy mode one gets an equation : 

"St + ^ = Sa, 
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where r\ is shear viscosity, S a is an acoustic source proportional to a gradient of acoustic intensity, index e 
concerns the entropy mode. An attenuation of acoustic wave leads to slow but continuous decrease of density 
so (^jf) = becomes incorrect as well as ^§f ^ = 0, since p = po + p a + p e . Therefore, if the acoustic heating 

caused by periodic acoustic wave is accounted, V£/ = becomes incorrect. 

These inconsistencies call new ways to theoretical investigation of acoustic streaming, including stream- 
ing caused by non-periodic acoustic wave. We plan to start with all conservation equations for compressible 
thermoviscous flow, then to define modes by establishing of the relations between perturbations (section 2), 
to demonstrate how the well-known linear and nonlinear evolution equations may be obtained by projecting 
(sections 3,4), and, finally, to apply projectors to calculating of the streaming caused by pulse acoustic wave 
(section 5). 

The idea to decompose the linear flow into specific modes is not novel and has been exploited for a long time 
, see [9] and referred papers, where homogeneous background with sources of heat and mass are considered. 
The concrete ideas of doing it automatically by using of projectors in the wide variety of problems (also flow 
over inhomogencous media like bubbly liquid [10], flows affected by external forces including the gravitational 
one which changes the background density and pressure [11], both one- and multi-dimensional problems, see 
[12] ). The principle advance is an expansion of the ideas into area of nonlinear flow: to get nonlinear coupled 
evolution equations for interacting modes and to solve the system approximately. The deriving of the final 
nonlinear system is algorithmic, one just acts by projectors at the initial system of nonlinear equations. Some 
ideas concerning separation and interaction of acoustic- gravity and Rossby waves may be found in [13]. 



2 Modes and projectors for the flow with equations of state in the 
general form. 

The mass, momentum and energy conservation equations for the thermoviscous flow read: 

| + V(ptf) = 



dv 



+ (vV)v 



= -Wp + r]Av+ (<r+ |) V(W) 



- + (vV)e 
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\dxk dxi 3 dxi 

Here, among already mentioned variables, p is pressure, e, T- internal energy per unit mass and temperature; 
<T, xare bulk viscosity and thermal conductivity coefficient (both supposed to be constants as well as ij ), Xi 
- space coordinates. Except of the dynamic equations (1), the two thermodynamic relations are necessary: 
e(p, p),T(p, p) . To treat a wide variety of substances, let us use the most general form of these relations as 
expansion in the Fourier series: 
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The background values are marked by zero, perturbations are primed, Cy means specific heat per unit mass 
at constant volume, E\, ...Oi, ...are dimensionless coefficients, 

(3) 



9i = i*£*, e 2 



where k - ± ( &) ,0 = ~ ± (U 

For the tree-dimensional flow v = (v x , v y ,v z ) in the Cartesian coordinates ~x = (x,y,z) over homogeneous 
stationary background without mean flow v = , the system (1), (2) looks: 

+ 1 - " Av ' - (" + 1) Yx^ v) = - po{ ™ )v * + itx - ^' Av * - ( ? + 1) p 'Tx^ v) 
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where Q,S are constants depending on the equations of state (2) : 

Q 

s 



±(-l + 2^—^E 3 + E 5 ), 



l-E 2 

For an ideal gas, Q = —7 = —Cp/C v , S = 0. 



1 1 - Eo 

(l + E 2 + 2E 4 + —-^E 5 ). 
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In order to simplify calculations, the geometry of beams will be considered. The equivalent system in the 
dimcnsionlcss variables : 



~v '*, a? \,p*,p*,U : ~v = ecv \,p' = ec 2 p p*,p' = sp p^,lr = (Xx*/</ji,\y*,\z*/y/jl),t = Af*/c, 



(6) 



(c = \J P " i p E^ '^ s ' adiabatic sound velocity, A means characteristic scale of the longitudinal perturbations, p 
is small parameter expressing the relation the longitudinal and transverse scales of perturbation, e is small 
amplitude parameter) looks as follows (asterisks for dimensionless variables will be later omitted): 



d_ 

dt 



V> + Lip = S(p + S(pt v , 



where ip is column of field perturbations 

IP = ( V X Vy V Z P p) 

L is a linear matrix operator : 
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dxdy 

^fpd/dx 
V Jpd/dx 

with dimensionless parameters 

si _ (C+a/3) r2 _ 77 rl _ xgi ^2 _ 

°1 pocA '"1 pocA'°2 pocAoC„Bi ' "2 pocA C„(l-B 2 ) ' 

Dimensionless operators V, A look : V = ( pd/dx d/dy pd/dz ) , A = pd 2 /dx 2 + d 2 /dy 2 + pd 2 /dz 2 , 
<p is the quadratic nonlinear column: 



X#2 



( -(v7)v x + yfppdp/dx \ 
-(v7)vy + pdp/dy 
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and iptv is the quadratic nonlinear column 0(/3) appearing in the viscous flow: 
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For a linear flow defined by the linearized version of system (7) 

d 



dt 



ip + = 0, 



(12) 



a solution may be found as a sum of planar waves: v x = v x ( k ) exp(iojt — ik ~x), ....where 

k = (k x , ky, k z ) 

is the wave vector. It is convenient to go to the Fourier transforms marked by tilde for the flows over the 
homogeneous background. In the Fourier space,— ik x means d/dx, iu means <9/<9i,and (6) yields in the five 
roots of dispersion relation ((3 = 5\ + 5\ + 8\ + <5f ): 



wi = fl + if3fl 2 /2, lu 2 = -ft + if3tt 2 /2, lu 3 = -i5 2 ft 2 , uj^ = i8 2 Q 2 , uj 5 = 2<5^ 2 ,where 



^1 Jvy ~\~ 
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(13) 



Also, an operator (—iky)" 1 represents in k -space operator J dy . Constant of integration should be chosen 
accordingly to the concrete physical problem. The first two roots relate to progressive (acoustic) modes of 
different directions of propagation, the third one relates to the entropy mode, and the fourth and fifth ones 
- to the vortical ones. For the real substances, (3 > 0,and, 8 2 < 0,that corresponds physically correct sign of 
imagine parts of all frequencies. The modes of linear flow are determined by relations of amplitudes of plane 
waves v x (k x ,ky 7 k z ) ,. . . Modes as eigenvectors of a linear problem in the Fourier space look: 
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All calculations of modes and projectors have accuracy up to the terms of order p, (3. As basic variable for 
the first three modes the perturbation of density is chosen, and for the last two- the stream function, since for 
the both vorticity modes density keeps unperturbed. So, there are five independent values indeed, one for every 
mode, that define other wave perturbations of the mode uniquely. Any field of the linear flow as a solution of 
a linearized equation (12) may be presented as a sum of independent modes. 

The next step is to get projectors that decompose a concrete mode from the overall field ip. Let us define a 
matrix M in the following way: 
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P3 
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ip = 'YjjJi = Mipi + M 2 p 2 + M 3 p 3 + M 4 </? 4 + M 5 tp 5 . 



(15) 



so that matrix M consists of five columns, 



M = ( Mi M 2 M 3 M 4 M b ). 

and the inverse matrix M^ 1 consists of five lines 

_1 = ( M- 1 M- 1 M- 1 M- 1 Mr 1 ) , 



(16) 
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where lines M 1 , ..M 5 1 are as follows: 
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(17) 



Accordingly to definition of matrix M 1 projectors may be determined, so as 

P^ = Vi,...,P 5 V = V'5. _ _ _ _ _ X 18 ) 

Relations (14), (15) written for the first mode as an example Mf tp — pi,ipi — M\p\ lead to P\ = Mi -M-f 
accordingly to (18) and so on for the all other modes: 



P i = M i -Mr\i = l,..5. (19) 

At any moment of evolution a concrete mode is distinguished from the overall field by the correspondent 
jjector. Pi 
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(20) 



Matrix projectors satisfy common properties of orthogonal projectors: 

J2 Pi = T, Pi-Pn =0iii^n, Pi-Pi = Pi if i = n, (21) 

i=i 

where 1,0 are unit and zero matrices. The inverse transformation of formulae (20) to the (~x,t) space may 
be easily undertaken. 
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3 Linear flow and evolution equations. 



A linear flow may be decomposed to modes uniquely accordingly to (18). Linear evolution equations for every 
mode may be originated from the linear version of (7) by acting of the corresponding projector: P ( JjV' + Lip) = 
0. Note that all projectors do commute with both d/dt I and L (I is unit matrix). When one acts by projector 
at the linear system (12), five equations for all the components of every mode appear. For example, an evolution 
equation for the dimensionless perturbations of density for both acoustic modes are as follows: 

%±^±f / A ±Ph2 dy-l %^ = 0, (22) 
where A ± = 1^ + 1^. 

An evolution equation for the axial velocity of the entropy mode looks: 

dv y3 X2 d 2 vy3 _ n (0 o\ 

Si 2Sy r ~ ~ U ' \ A6 > 

and for the axial velocity of two vortical modes is: 

(24) 

These linear evolution equations arc well-known .Equations (21)-(23) are explicit up to terms of order 
0(f3y/JI, p 2 ) due to accuracy of evaluated projectors. In the paper [9], equations (22)- (24) are presented as three 
equations for the three basic modes, since there was no subdivision of acoustic and vortical modes into two 
independent branches everyone. Also, the evolution equations from [9] involves sources of mass, momentum 
and heat that easily may be added in to the right-hand side of the initial system (7) as well as to the equations 
(22)-(24). 



4 Nonlinear flow and coupled evolution equations. 

At this point, we declare that relations inside every mode become fixed in weakly nonlinear flow. Acting by 
projectors at the original system (7) with non-zero nonlinear part leads to coupled equations for interacting 
modes: P ( Jr^ + Lip) = eP (tpi + (fitv)- Modes become separated in the linear side, and one should account 
that the overall field is a sum of all modes inputs in the nonlinear side. The proper returning to the ~x 
space should be proceeded as well. To illustrate the possibilities of the projecting method, some illustrations 
concerning the famous evolution equations deriving are presented below. 

4.1 One-dimensional nonlinear flow. 

Let k x = 0, k z = 0, that corresponds to one-dimensional flow along y-axis. For the planar geometry, only three 
modes exist: two acoustic and the entropy one. Acting by Pi at the system (5), one get an evolution equation 
for the rightward acoustic mode consisting of three equations for , pi , pi , and v\ . An evolution equation for the 
non-dimensional density perturbations of the rightward acoustic mode is as follows: 



dpi dpi (3 d 2 V! = e ( v d v { ^_ _ V JL + \n + s } —1 
dt dy 2 dy 2 2 \ dy dy dy dy 

Here, p, p, and v in the right-hand part represent a sum of all modes: 

V = Vi + V 2 + V 3 , p = Pl+p 2 +P3, 

P = Pl+P2 + P3- (26) 

All cross nonlinear- viscous terms are neglected, as well as higher order nonlinear ones. If the nonlinear terms 
of the first mode are kept only, we get the known Burgers evolution equation, going to the famous Earnshow 
one in the limit /3 — > : 



£ n Pi—. Pi - =0- (27) 



dpi dpi -Q-S+l d (3d 2 Pl 
~dT + ~dy~ + £ 2 Pl dy~ Pl ~2~dy~ 2 

Physically, this case relates to a dominant rightward acoustic mode. For an ideal gas, again, Q = —7 = 
—Cp/C v , 

s = o. 
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4.2 Dynamics of beams. 



In the same way, the famous Khokhlov-Zabolotskaya equation [14,15] (KZ, in the limit j3 — > 0) and Khokhlov- 
Zabolotskaya-Kuznetsov equation (KZK) for rightward (leftward) propagating beams follows when acting by 
P\P2 at the system (7): 

dpi,2 , dp 1>2 , p f . , /3d 2 p 1 . 2 , -Q-S + l d 

~dT ~dy~ 2 J A±Pl ' 2dV -~2^F ± £ 2 P ^dy PX * = °' (28) 

This evolution equation may be understood as a limit for the progressive mode self-action, meaning that only 
this mode nonlinear terms are kept. It seems be suitable definition though the 'self-action' is usually reserved 
for the phenomenon of acoustic heating when acoustic mode causes the entropy mode during its propagation 
thus changing the surrounding [8] . 



5 ACOUSTIC STREAMING CAUSED BY IMPULSES. 



There are some reasons to find another approach to acoustic streaming. First, the procedure of temporal 
averaging is a basis of the modern theory of acoustic streaming and all corresponding results relate to averaged 
fields. All non-periodic acoustic sources thus are left of account. Second, the problem of modes interaction looks 
more extended. Actually there exist three types of independent modes. The initial amplitudes of these modes 
are determined by initial conditions of a concrete problems and may be evaluated by projecting of the overall 
initial field into every modal field. The acoustic streaming imposes an ability of predominate acoustic mode. The 
acoustic field causes slow varying with time flow which may grow with time since nonlinear effects can store, that 
is namely the acoustic streaming. So, the acoustic streaming relates to one type (though important) of possible 
nonlinear interactions in flow, namely to the beginning of acoustic mode evolution. Later, an amplitude of the 
vortical mode (and the entropy one also) may become so large that one can not treat these modes as secondary 
and other approach for nonlinear interaction of modes should be developed. Let us mark, that even traditional 
subdivision of overall flow into 'slow' and 'quick' components eliminates the entropy mode unclcarly. It is also 
'slow' but essentially possesses non-zero density perturbation which slowly varies with time ( dp 3 /dt = in 
linear approach), so the procedure of temporal averaging of continuity equation fails with storing of nonlinear 
effects. Also, there are many more complicated flows over inhomogeneous media and/or with background flows 
that may be algorithmically solved by the pointed method, see conclusion. 

We will show, at first, that projecting results in the well-known equations for an acoustic streaming in terms 
of V , and, the second and more important, that there appear the new equations for acoustic streaming caused 
by impulse (non-periodic) acoustic signals. 

To consider interactions of vortical and acoustic modes, one acts by projector P4 + P5 at the system (7) 
to get an evolution equation for vorticity ip^ + 935 . The right-hand nonlinear vector involves, again, all modes 
inputs. Then, to calculate the nonlinear interactions of only vortical and rightward propagating acoustic mode, 
the right-hand vector should be thought as a sum of vortical and the first acoustic mode: 



l-BA ± fdyfdy-l% 

V 1 J 



Pi 



( V X 4 \ 




V / 



/ \ 

Vy 5 

v z5 


V / 



(29) 



An operator Aj_ J dy J dy represents in ~x -space operator k *£i kz .Acting by P 4 + P 5 at the system (5) leads 

V 

to an evolution equation (27) for vorticity ip(x ,t), which we rewrite going to velocity components V x — v X 4 — 
-dipi/dy, V y = v y4 + v y5 = ^/Jidtpi/dx + y/JId(p 5 /dz, V z = v z5 = -d<p 5 /dy : 



dV x 2 d 2 V x 
dt 1 dy 2 



-e(VV)V x +eF lx 



(30) 



where in right-hand side of (30) only two term are left: the first relating to the vortical mode, and the second 
F\ consisting of quadratic inputs of the first progressive acoustic mode: 
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F i* = (l-»&JdyJdy -y/Ji&fdy -p^jdyjdy ) (<pi + <pitv), 



(31) 



where <pi,tpit v are supposed to consist of specific perturbations for the rightward dominant acoustic mode, 
see (10). For example, <pi looks 



<Pi 



\ 



There are all other possible cross nonlinear terms in the right- 



-(viV)vi x + yfjipxdpi/dx 
-(viV)v-iy + pidpi/dy 
-(viW)vi z + yfjlpidpi/dz 
[Qpi + S,oi](Wi)-_(i/itf)pi 
-pi(Wi) - {viV)pi 

hand side of (30) which are out of interest in the present investigation since the rightward acoustic mode 
supposed to be dominant at least at the initial stage of the evolution. Writing on the terms of order not higher 
then 0fj, , one goes to the expected result: 



V 



Jdyjdy ) (pi = 0, 



that means that acoustical streaming may exist only in the thermoviscous flow. Finally, F\ x is as follows 
F lx =(l-p£ 7 JdyJdy -y/Ji&fdy -^Jdyjdy ) <p Uv = 



(32) 



V7* 



( dpi dpi 
\ dx dy 



dy Pi 



d 2 Pl 
dy 2 



(33) 



To show that both effects play role: (1) presence of the thermoviscous terms in tpi, and (2) presence of 
thermoviscous terms in ipu v , corresponding terms are included in square brackets, first and second ones though 
they partially compensate each other. Finally, (33) goes to the following: 



/ dpi dpi 



0.5p : 



d 2 Pl 

dxdy 



(34) 



and the constant of integration of the last term should be chosen accordingly to the real initial perturbations 
and /or boundary regime. It is remarkable that only three first elements of ipitv participate in the final radiation 
force, namely that three ones are responsible for the losses of momentum in the flow. The forth element expresses 
the thermoviscous loss of energy. This term (and also the fifth one) plays role in the forming of the heating, 
see the last line of P 3 given by (20). Relations of the specific perturbations of eigenvector ipi (14) are used to 
express F\ x through p\. 

The transverse component of the force caused by sound F\ x given by formula (34) does not possesses terms 
depending on the other transverse coordinate z because these terms are of order 0p?l 2 . In all calculations, the 
properties of specific modes and related projectors applied only. The temporal averaging was not proceeded. It 
has been proved that the acoustic force should be rotational [16]. In the referred papers, the rotational part of 
the averaged is just selected from the general expression. Projectors due to their properties yield in the 

vortical forces automatically. 



5.1 Limit of the quasi-periodic acoustic source 

To compare formula (30) with this caused by periodic at the tranducer source, let us take the transverse acoustic 
force (radiation force) given by Gusev, Rudenko [17] 

= VMi|^ 2 /2) exp(-/3y) (35) 

rewritten in our variables; note also that longitudinal coordinate there is x instead of y. The radiation force 
(35) relates to the streaming evolving in the sound field of a highly attenuated beam with a quasi-periodic 
source as follows: 
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Pi (x, V,t) = p a O(x) exp(- ^y) sin(i - y) 
see [12] for more details. Calculating based on formula (34) yields in the following: 



1 + - sin(2(i • 
4 



v)) + o((3 2 ) 



(36) 



(37) 



Averaging over integral number of period of the sound wave in the leading order gives {F\ x ) = &i x . In 
contrast to the known results, (34) fits to calculate force caused by any acoustic source including non-periodic 
ones. Moreover, evolution of streaming given by temporal averaging is not suitable in principle for detail tracing 
with small temporal step because of the very initial procedure of averaging over time whose interval is much 
longer than a period of a sound wave. 



5.2 Radiation force cased by mono-polar acoustic source 

As an appropriate acoustic source satisfying the limitations of [17], let us take a mono-polar two-dimensional 
wave as follows [18]: 



Pi{x,y,t) 



2/3 . 2 , 
exp(— x ) — 



exp i 



720 



(38) 



UP [C - Erf(r/V^) 
®~ 2 S+1 1S a parameter of nonlinearity, see (5). For ideal gases, e = Constant 



where £ = (iy, t = t — y, e 

C is responsible for the shape of the curve: large C provides a curve close to the Gauss one. The function 
Fi x achieves negative minimum at the transversal distance x m = y/2/2 and is an odd function as follows from 
(34), (38). Results of the calculating of in accordance to formula (34) arc shown at the figures 1,2 at the 

transversal point x = x m . Figure 1 shows the dependence of ^j=^ on y at t = 1, and figure 2 shows the same at 
t = 3. Constants of mono-polar source (35) used in calculations are: C = 2, e = 1.2, (3 = 0.1. Also, the constant 
of numerical integration of the last term of (31) is chosen to the radiation force be zero at y— > oo. 

All pictures show that streaming develop with some delay after the source passing. The reason is non-local 
relation between the source and induced streaming. The velocity field may be calculated accordingly to the 
formula (30) completed with equations for the other components of V. At the first stage of evolution, the 
nonlinear term ( V V) V may be neglected and the linear equation like this of thermal conductivity with viscous 
coefficient cased by shear viscosity df should be solved. The acoustic source in the right-hand side of (30) is 
namely the radiation force. The temporal behavior of the radiation force at the two longitudinal coordinates 
y = 1 and y = 3 is presented at the figures 3,4, x — V2/2. The mono-polar source is naturally attenuated during 
its propagation along y-axis. Note that the radiation force is multiplied by the large value -^=g at all pictures, 
so the real scale is much less. 




Fig.l Longitudinal distribution of the radiation force ^^(y) (thin line) and the acoustic source pi(y) (bold 

line) at t = l,x = V2/2. 



9 




Fig. 2 Longitudinal distribution of the radiation force 4=^(2/) (thin line) and the acoustic source pi(y) (bold 

line) at t = 3,x = y/2/2. 




Fig. 3 The temporal distribution of the radiation force -^^(t) (thin line) and the acoustic source pi(t) (bold 

line) at y = l,x = y/2/2. 




Fig. 4 The temporal distribution of the radiation force -jj^(t) (thin line) and the acoustic source p\{t) (bold 

line) at y = 3, x = V2/2. 

6 Conclusions 



The basic idea of the paper is to separate modes accordingly to their properties in the weakly nonlinear flow. 
At first, modes as eigenvectors of the linear flow should be defined. In the other words, the relations of 
the specific perturbations inside every mode should be established. Both homogeneous and inhomogeneous 
backgrounds ( see paper on interacting modes in bubbly liquid [10]), media affected by external forces including 
the gravitational one which changes the background density and pressure [11] are treated in the algorithmic 
way. The definition of modes is unique, determined by the linearized differential conservation equations only. 
Proper initial or/and boundary conditions are given by the proper superposition of modes. 

The second step is to get coupled nonlinear equations for the interacting modes and to solve it approximately. 
If one of the acoustic modes is dominant the generation of the entropy mode and the other acoustic mode in 
the plane geometry (namely heating and the reflected wave) is governed by the simple analytical formulae. The 
evolution equations may be corrected up to the higher order nonlinear terms due to increasing influence of the 
other generated modes [10,12]. 

The principal advance in comparison to the widely used approaches consists in possibility to treat any 
sources: initial mixture of all possible types of motion, non-periodic acoustic source and so on. In the present 
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paper, radiation force caused streaming is calculated for any acoustic source, an example of the monopolar pulse 
source is considered, the limit of periodic source is traced. 
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